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ABSTRACT 


During  the  months  of  July  through  September,  1968,  a 
series  of  atmospheric  nuclear  explosions  was  detonated  by 
France  in  the  South  Pacific.  The  events  analyzed  are  those 
of  July  7,  July  15,  August  3,  August  24,  and  September  3, 
1968,  The  analysis  includes  matched  filtering  for  relative 
amplitude  (H)  measurements,  body  and  surface  wave  magnitude 
comparisons  phased  summing,  and  ARZ-ERZ  measurements.  The 
events  are  discussed  individually,  and  a  summary  of  per¬ 
tinent  data  is  presented  in  several  figures  anu  tables. 
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INTRODUCTION 


During  the  months  of  July  through  September,  1968  a 
series  of  atmospheric  nuclear  explosions  was  detonated  by 
France  in  the  South  Pacific.  The  events  from  that  series 
unaxyzed  in  the  report  are  those  of  July  7,  July  15, 

August  3,  August  24,  and  September  8,  1968.  The  analysis 
includes  matched  filtering  (Alexander  and  Rabenstine,  1967a, 
1967b)  for  relative  amplitude  4S>  measurements,  body  and 
surface  wave  magnitude  comparisons  (Lambert  and  Turnbull, 
1968),  phased  summing,  and  ARZ-ERZ  measurements  (Turnbull 
and  Lambert,  1968).  The  events  will  be  discussed  individually 
m  the  following  sections,  with  a  summary  of  pertinent  data 
presented  in  several  figures  and  tables.  (  V 


Data  Anali 


The  five  events  of  interest  occurred  essentially  at 
the  same  coordinates.  The  data  we  wished  to  analyze  were 
those  recorded  at  LASA,  TFO ,  UBO,  and  five  LRSM  stations 
(LCNM,  MNNV,  PGBC,  RKON,  and  SV3QB) .  The  data  from  TFO 

and  UBO  proved  unusable  because  neither  LP  array  was  yet 
operational. 

The  LASA  data  were  recorded  digitally  at  5  samples 
per  second  and  decimated  to  1  sample  per  second  for  analysis 
The  LRSM  data  were  recorded  on  analog  tape,  digitized  at 
the  SDL  at  5  samples  per  second,  and  then  decimated  to  1 
sample  per  second  for  analysis. 

The  LASA  data  were  demagnified  using  results  of  .04  Hz 
sine  wave  electro-magnetic  calibrations  recorded  in  March  1967 


Obviously  these  calibrations  are  not  accurate  for  current 
data,  but  they  are  the  most  recent  that  we  have  available. 
This,  however,  should  not  affect  the  relative  amplitude 
(ct)  values.  The  LRSM  data  were  demagnified  by  comparison 
with  the  film  records  and  indicated  magnifications. 

Only  the  sensors  in  the  A,D,E,  and  F  subarrays  at  LASA 
were  used,  in  order  to  assure  incoherent  noise  (Alexander 
and  Rabenstine,  1967b). 

The  primary  technique  of  our  analysis  was  the  matched 
filter.  A  search  for  events  which  originated  in  the  same 
source  region  with  sufficient  magnitude  to  generate  a  good 
dispersed  wavetrain  was  undertaken  using  the  library  pro¬ 
gram  DEPTHMAG.  The  source  region  was  shown  to  be  extremely 
aseismic,  with  only  one  event  meeting  the  needed  specifioa- 
tions  (i.e.,  within  10°).  An  earthquake  occurred  at  18. 4S, 
132. 9W  of  magnitude  5.5  on  March  6,  1965.  LASA  was  not  in 
operation  at  that  time.  However,  an  ideal  matched  filter 
was  furnished  by  one  of  the  explosions  itself;  the  event 
of  August  24.  The  body  wave  magnitude  was  4.8  and  the 
surface  wave  magnitude  was  found  to  be  4.7.  A  dispersed 
wavetrain  with  good  S/N  was  readily  visible  on  the  records 
at  LASA.  Thus,  the  a  values  measured  are  amplitudes  of  the 
various  events  relative  to  the  event  of  August  24. 

The  bandpass  filter  which  was  used  on  the  raw  data 
either  before  phased  summing  or  when  matched  filtering  was 
simply  a  shaping  applied  to  the  Fourier  transform.  The 
transform  was  multiplied  by  a  linear  taper  which  went  from 
1.0  to  0.0  as  the  frequency  went  from  0.02  to  0.0167  Hz  and 
0.1  to  0.125  Hz.  Transform  values  outside  the  frequency 
range  of  0.0167  to  0.125  Hz  were  set  to  zero. 
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A.  July  7 


Long  period  analysis  of  this  event  by  all  techniques 
proved  quite  futile.  Approximately  one-half  hour  before 
detonation,  an  earthquake  of  magnitude  5.1*  (m^)  occurred 
in  Mindanao,  Phillipines  (9.6N,  126. 5E).  The  paths  were 
such  that  the  wavetrains  of  the  earthquake  and  the  explo¬ 
sion  overlapped  at  LASA  and  most  LRSM  stations.  Results 

of  the  phased  sum  are  shown  in  Figure  1  and  of  the  matched 
filter  and  phased  sum  in  Figure  2. 

b.  July  is 


This  event  was  determined  to  have  a  body  wave  magni¬ 
tude  (mb)  of  4.3*  and  a  surface  wave  magnitude  (hs)  of  4.1. 
The  surface  wave  magnitude  is  considerably  larger  than  that 
for  an  underground  nuclear  explosion  of  comparable- body 
wave  magnitude  (approximately  two  orders  of  magnitude),  an 
indication  that  an  atmospheric  nuclear  explosion  -generates 
a  great  deal  more  long  period  energy.  In  a  plot  of  Ms  vs. 
mb  (Figure  13),  the  events  of  this  series  are  shown  along 
with  earthquakes  and  underground  nuclear  explosions.  For 
the  same  order  body  wave  magnitude,  we  find  that  atmospheric 
explosions  generate  even  mors  long  period  energy  than  earth¬ 
quakes.  The  energy  difference  is  further  displayed  in  plots 
of  ARZ  vs.  mb  (Figure  14)  and  ERZ  vs.  mb  (Figure  15). 

The  signal  as  recorded  at  LASA  was  readily  visible 
(Figure  3).  The  records  of  eleven  long-period  instruments 


All  body  wave  magnitudes  are  from  the  LASA  bulletin;  all 
surface  wave  magnitudes  are  computed  in  the  ARZ-ERZ  program 
(Turnbull  and  Lambert,  1968). 
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were  band-pass  filtered  from  10  to  50  seconds  (period)  and 
phase  summed  using  a  velocity  of  3.75  km/sec-  From  an  average 

S/N  ratio  of  4.0  for  each  channel,  the  sum  channel  S/N  was 

9.5. 


The  LASA  matched  filter  output  was  quite  impressive, 
as  could  be  expected  from  the  visible  signal  (Figure  4). 
ihe  S/N  ratio  of  ne  matched  filter  output  had  an  average 
of  11.2  (a  avg.  =  0.248),  and  the  S/N  ratio  of  the  sum  of  the 
matched  filter  outputs  from  eleven  channels  was  34.7 
(a  -  0.238).  The  summed  matched  filter  output  was  quite 
symmetric  and  narrow  indicating  high  correlation  between  the 

filter  and  event.  The  standard  deviation  of  the  S  values  was 
14%  of  the  mean  value. 

When  the  five  LRSM  signals  were  matched  filtered  and 
beam  formed  with  the  LASA  matched  filter  outputs  on  4  km/sec, 
a  further  S/N  enhancement  to  47.1  was  achieved  (Figure  5), 
with  an  £  of  0.209  on  the  sum.  The  mean  of  all  S.  values  was 
0.234  with  a  standard  deviation  of  16%. 

C.  August  3 

This  event  was  determined  to  have  a  body  wave  magnitude 
(mfc)  of  4  1  and  a  surface  wave  magnitude  (Ms)  of  4.3.  As  in 
the  July  15  event,  a  great  deal  more  long  period  energy  than 
in  an  underground  explosion  of  the  same  m^  is  generated.  How¬ 
ever,  the  signal  as  recorded  at  LASA  was  not  as  discernible 
as  the  July  15  event  (Figure  6),  the  average  S/N  being  3.1 

and  the  S/N  of  the  phase  summed  trace  of  eleven  channels  only 
3.8. 

This  reason  for  the  very  small  indicated  increase  in  S/N 
is  probably  that  the  signal  levels  measured  on  the  individual 
traces  were  too  high  because  of  noise  being  mistakenly 
measured  as  signal. 
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The  LASA  matched  filter  output,  though,  showed  a  greater 
improvement  upon  summation  (Figure  7).  While  the  average 
S/N  for  eleven  matched  filtered  channels  was  still  only  3.1 
(S  avg.  =  0.0695),  the  sum  S/N  was  8.3  (a  =  0.06).  The 
standard  deviation  of  a  values  was  32%,  again  due  to  the 
relatively  low  S/N  ratios. 

Once  again  beamforming  the  5  LRSM  matched  filtered 
signals  with  the  LASA  matched  filtered  signals  on  a 
velocity  of  4  km/sec  increased  the  S/N  to  9.56  (Figure  8) 
with  a  sum  trace  a  of  0.057.  The  mean  of  all  £  values 
was  0,0698  with  a  standard  deviation  of  24%. 

D .  August  24 

This  event  had  a  body  wave  magnitude  (m^)  of  4.8 
and  a  surface  wave  magnitude  (Mg)  of  4.75.  The  Rayleigh- 
wave  train  was  very  well  defined  on  the  LASA  records, 
hence  its  use  as  the  matched  filter  for  the  other  events 
(Figure  9).  The  S/N  for  an  average  band  pass  filtered 
channel  was  17.9,  while  the  S/N  of  the  phased  sum  of 
thirteen  channels  is  51.7. 

E .  September  8 

This  event  was  the  second  largest  of  the  series  and 
had  a  body  wave  magnitude  of  4.6.  The  surface  wave  magnitude 
was  measured  to  be  4.6,  and  the  signal  at  LASA  was  well 
recorded.  The  previous  conclusions  on  surface  wave  energy 
from  atmospheric  explosions  are  again  verified. 

The  beam  forming  of  the  band  pass  filtered  surface  waves 
at  LASA  worked  quite  well  giving  an  improvement  in  S/N  from 

a  mean  of  6.3  on  the  individual  traces  to  14.1  on  the  sum 
trace  (Figure  10.). 
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The  matched  filter  also  worked  well,  as  expected  with 
a  large  signal.  The  LASA  matched  filtered  traces  achieved 
a  mean  S/N  of  13.51  (a  avg.  =  0.541)  and  summed  to  give  a 
S/N  of  37.77  (a  =  0.508)  (Figure  11).  The  LRSM  matched 
filtered  data,  which  had  mean  S/N  of  11.35  (a  avg.  =  0.449), 
combined  with  the  LASA  data  to  give  a  sum  S/N  of  42.33 
(a  =0.488)  (Figure  12).  The  standard  deviation  of  all  the 
values  was  22%  of  the  mean. 
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CONCLUSIONS 

1.  In  comparing  atmospheric  explosions  with  underground 
explosions,  we  find  the  former  generate  more  long  period  energy 
than  the  latter,  with  the  surface  wave  magnitudes  generally 
two  orders  of  magnitude  greater. 

2.  The  "best"  matched  filter  proved  to  be  a  large 
atmospheric  shot  from  the  same  region. 

3.  The  LRSM  matched  filtered  data  was  beam  formed  with 
the  LASA  matched  filtered  data  to  give  S/N  values  very  near 
that  which  is  predicted  for  coherent  signals  in  incoherent 
noise.  The  exception  was  a  case  where  the  actual  values  of 
S/N  on  the  individual  LASA  channels  were  too  low  to  measure 
accurately . 

4.  Reasonably  consistent  values  of  a  were  obtained  across 
the  network  of  stations  used. 
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Figure  4.  Event  of  15  July  1968.  LASA  LPZ  Rayleigh  waves 
'patched  filtered  with  signals  from  event  of 
24  August  1968.  Sum  of  LASA  matched  filtered 
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Figure  5.  Event  of  15  July  1963.  LRSM  LPZ  Rayleigh  waves 
matched  filtered  with  signals  from  event  of 
24  August  1968.  Sum  of  LASA  and  LRSM  matched 
filtered  signals. 
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Figure  9.  Event  of  24  August  1968.  LASA  LPZ  Rayleigh  waves 
band  pass  filtered  (10—50  sec.  period)  and  beam- 
formed  on  3.75  km/sec. 


Figure  10.  Event  of  8  September  1968.  LASA  LPZ  Rayleigh  waves 
band  pass  filtered  (10-50  sec.  period)  and  beam- 
formed  on  3.75  km/sec. 
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Figure  12.  Event  of  8  September  1968.  LRSM  LPZ  Rayleigh  waves 
matched  filtered  with  signals  from  event  of  24 

August  1968.  Sun  of  LASA  and  LRSM  matched  filtered 
signals. 
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Figure  14.  ARZ  vs  comparison  of  various  seismic  signals 


Figure  15.  ERZ  vs  xr^  comparison  of  various  seismic  signals 
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